Astronomy & Astrophysics manuscript no. virus-preprint 


©ESO 2012 


September 4, 2012 





Mapping Luminous Blue Compact Galaxies with VIRUS-P: 

morphology, line ratios and kinematics 

L. M. Cair6s\ N. Caon^'^ B. Garcia Lorenzo^'^ A. Kelz\ M. Roth^ P. Papaderos^, and O. Streicher^ 

^ Leibniz-Institut fiir Astrophysik, An der Sternwarte 16, D- 14482 Potsdam, Germany 

e-mail: luzma; akelz; mmroth; ole@aip.de 
^ Instituto de Astrofisica de Canarias, E-38200 La Laguna, Tenerife, Spain 
^ Departamento de Astrofisica, Universidad de la Laguna, E-38206 La Laguna, Tenerife, Spain 

e-mail: nicola . caon ; bgarcia@iac . es 
^ Centro de Astrofisica and Faculdade de Ciencias, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 

e-mail: papaderos@astro . up . pt 

Received May 15, 2012; accepted August 29, 2012 

ABSTRACT 

Context. Blue Compact Galaxies (BCG) are narrow emission-line systems, undergoing a violent burst of star formation. They are 
compact, low-luminosity galaxies, with blue colors and low chemical abundances, that offer us a unique opportunity to investigate 
collective star formation and its effects on galaxy evolution in a relatively simple, dynamically unperturbed environment. Spatially 
resolved spectrophotometric studies of BCGs are essential for a better understanding of the role of starburst-driven feedback processes 
on the kinematical and chemical evolution of low-mass galaxies near and far. 

Aims. We carry out an integral field spectroscopy (IFS) study of a sample of luminous BCGs, with the aim to probe the morphology, 
kinematics, dust extinction and excitation mechanisms of their warm interstellar medium (ISM). 

Methods. IFS data for five luminous BCGs were obtained using VIRUS-P, the prototype instrument for the Visible Integral Field 
Replicable Unit Spectrograph, attached to the 2.7 m Harlan J. Smith Telescope at the McDonald Observatory. VIRUS-P consists of a 
square array of 247 optical fibers, which covers a 109'' x 109'' field of view, with a spatial sampling of 4'.' 2 and a 0.3 filling factor. 
We observed in the 3550-5850 A spectral range, with a resolution of 5 A FWHM. From these data we built two-dimensional maps 
of the continuum and the most prominent emission-lines ([On] /13727, Hy, H/? and [Oiii] /15007), and investigate the morphology 
of diagnostic emission-line ratios and the extinction patterns in the ISM as well as stellar and gas kinematics. Additionally, from 
integrated spectra we infer total line fluxes and luminosity- weighted extinction coefficients and gas-phase metallicities. 
Results. All galaxies exhibit an overall regular morphology in the stellar continuum, while their warm ISM morphology is more 
complex: in II Zw 33 and Mrk 314, the star-forming regions are aligned along a chain- structure; Haro 1, NGC 4670 and III Zw 102 
display several salient features, such as extended gaseous filaments and bubbles. A significant intrinsic absorption by dust is present 
in all galaxies, the most extreme case being III Zw 102. Our data reveal a manifold of kinematical patterns, from overall regular gas 
and stellar rotation to complex velocity fields produced by structurally and kinematically distinct components. 
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1. Introduction 

Blue Compact Galaxies (BCGs) are low-luminosity (Mb ranges 
from -12 to -21) and low-metallicity (7.0 < 12 -h log(0/H) < 
8.4) systems, which have optical spectra similar to those of H ii 



region s in spiral galaxies (Thuan & Martin 1 981t|Kunth & Qstlin 
[2000 ). With star formation rates between 0.1 and 1 Mq yr~^ 
( [Fanelli et al.|1988[ [Hunter & Elmegreen|2004l ), they represent 
the most extreme manifestation of star-forming (SF) activity in 
late-type dwarf galaxies in the local Universe. 

BCGs are of paramount importance in extragalactic astron- 
omy and observational cosmology research, as they hold key 
clues to understand fundamental topics such as galaxy formation 
and evolution, and the ignition and self-regulation of SF activity 
in relatively isolated late-type galaxies. 

They represent an unparalleled link to the early Universe, 
since in the well accepted framework of a Cold Dark Matter 
(CDM) Universe, structure formed hierarchically, with small 
scale objects (low mass halos) collapsing first at relatively large 
redshifts. About two percent of the BCGs have very low metal- 



licities, (12 -\- log(0/H) < 7.65); these objects, often referred 
to as extre mely Metal Deficient (XBCD) galaxies (Papaderos] 
et al.|2Q08] ), constitute the best local analogs to the distant sub- 



galactic units from which larger galaxies are formed (Papaderos 



|etaL| 1998 ; Kniazev et al.|2004 ). Also, the more luminous BCGs 
(or different subsamples of them) have been regarded as the lo- 
cal counterparts of different high redshift SF galaxy populations 
d Mas-Hesse et al.||2003[ [Grimes et al.||2Q09| ). Detailed analysis 
of these nearby objects are essential to interpret the observations 
of the more distant SF galaxy population, as their proximity al- 
lows for studies focused on their stellar content, kinematics and 
abundances, with an accuracy a nd spatial resolution that ca nnot 
be achieved at higher redshifts ( Cairos et al.|2QQ9a]b 2010). 

Moreover, low-luminosity systems are the best laboratories 
to study the star- formation process. Without spiral density waves 
to act as a triggering mechanism for star formation, these galax- 
ies enable us to investigate in an unbiased manner all of the other 
factors that trigger and control SF activities. In addition, BCGs 
offer a unique opportunity to test star-formation models at very 
low metallicities, similar to those of the early Universe. 
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In spite of the strong interest gained for BCGs in the last 
decades, several key issues in the field, namely, the evolutionary 
status of BCGs, their Star Forming History (SFH) or the mecha- 
nisms that trigger their star-forming activity, are still open. 

One of the reasons is probably that most of the BCG work 
published so far deal with analysis of samples by means of sur- 
face photometry ("Papade ros et al.||1996 |Cair6s et al.||2001a|b[ 



Doublier et aDll997, 19991 |Gil de Paz et al.||2QQ3t |Gil de Paz| 



& Madore 2005 ), and spectroscopic studies typically only cover 
the brightest SF regions of these systems. Studies combining sur- 
face photometry with spatially resolved spectroscopy and spec- 
tral synthesis, and aiming at a quantitative study of the formation 
history of stellar populations in BCGs, are scarce. 

That very few spectrophotometric analyses can be found in 
the literatur e, and vir tually all of them focused on one single 
object (Guseva et al72003a b cl [C airos et al.|2002l[2007]), is es- 
sentially due to the large amount of observing time that con- 
ventional observational techniques require: acquiring images in 
several broad-band and narrow-band filters, plus a sequence of 
long- slit spectra sweeping the region of interest, translate into 
observing times of a few nights per galaxy, making analysis of 
statistically meaningful samples not feasible in terms of observ- 
ing time. Moreover, such observations usually suff'er from vary- 
ing instrumental and atmospheric conditions, which make com- 
bining all these data together complicated. 

Integral Field Spectroscopy (IFS) off'ers a way to approach 
spectrophotometric BCG studies in a highly effective manner 



(Izotov et al. 2006a[ |Garcia-Lorenzo et al.| 


2008 IVanzi et al.| 


2008[ 


James et al.||2009 


Cairos et al.||2009a 


b 2010]). IFS pro- 



vides simultaneous spectra of each spatial resolution element, 
under identical instrumental and atmospheric conditions, which 
is not only a more efficient way of observing, but also guaran- 
tees the homogeneity of the dataset. From only one Integral Field 
Unit (IFU) frame we can produce a series of broad-band images 
in a large number of filters, narrow-band images in a large set of 
bands, as well as derive a large collection of observables from 
the spectrum (e.g. line indices). In terms of observing time, IFS 
observations of BCGs are one order of magnitude more efficient 
than traditional observing techniques, implying that now spec- 
trophotometric studies of large samples of BCG galaxies have 
become viable. 

The work presented here is part of the IFS -BCG survey, a 
challenging long-term project that aims to perform an exhaustive 
spectrophotometric survey of a large sample of BCGs by means 
of IFS. The main scientific goals of this projects include: i) to 
disentangle and characterize the diff'erent stellar populations in 
BCGs (e.g. constrain their ages. Initial Mass Function and metal- 
licity); ii) to investigate the evolutionary status of the galaxies 
and constrain their SFH; iii) to probe the mechanism responsible 
for the actual burst of star formation; iv) to analyze the feedback 
eff'ects between massive stars and the Interstellar Medium (ISM) 
in dwarf galaxies; v) to investigate the recent suggestion that 
some BCGs harbor Active Galactic Nuclei (AGN) ass ociated 
with intermediate mass black holes [Izotov et al.| ( |2010| ); vi) to 
provide an accurate dataset of photometric and spectroscopic pa- 
rameters for a large sample of nearby SF dwarf galaxies, the es- 
sential template to understand the results of the investigations at 
intermediate and high-z. 

This is the fourth in a series of papers presenting the find- 
ings of the project, in which we report on results derived for five 
luminous BCGs observed with VIRUS-P. In the first and sec- 
ond papers of the series ( [Cairos et al.|2009a|b ) we illustrated the 
full potential of this study by showing results on two represen- 
tative BCGs, Mrk 1418 and Mrk 409, both observed with the 



Potsdam multi- aperture spectrophotometer (PMAS), attached at 
the 3.5-m telescope at Calar Alto Observatory. A more detailed 
description of the project objectives, IFU observation and reduc- 
tion techniques and results f or another eight obj ects observed 
with PMAS, was presented in [Cair6s et al.[ ( [2010l ). 

This paper is structured as follows: In Sect. |2] we describe 
the observations, the data reduction process and the method em- 
ployed to derive the maps. In Sect. |3] we present the main out- 
comes of the work, that is, the flux, line ratio and velocity maps, 
as well as the results derived from the analysis of the integrated 
spectra. Finally, the main findings for each galaxy are discussed 
in Sect. [4] and summarized in Sect. jS] 



2. Observations and Data Processing 

2.1. Galaxy sample and observations 

The whole survey includes the mapping of about 40 galaxies, 
chosen as to cover the wide range of luminosities (-12 > Mb > 
-21) and metallicities (7.0 < 12-Flog(0/H) < 8.4) found among 
BCDs. We took also care that the diff'erent morphological classes 
observed among BCGs be well represented in the sample ( Cairos] 
|etalJ2001a|). 

In this paper we present results for fiv e luminous BCGs 
( jBergvaU & Qstlin|2002[ [Cairos et al.|2002| ). These objects are 
the most luminous (Mb < -17), and typically also the more ir- 
regular BCGs; they tend to have an important contribution of 
older stars, significa nt amounts of dust and very distor ted star- 
burst morphologies ( Cairos et al.[[2003t Bergvall[[2012[ ). As we 
have already stated in the introduction, these luminous BCGs 
have a major relevance in cosmological studies, as they (or dif- 
ferent subsamples of them) have been regarded as the local coun- 
terparts of diff'erent galaxies populations at higher redshifts (z ~ 
2-6), as Lyman-break g alaxies ( Grimes et al.|2009] ) or Lyman-a 
emitters ( [Mas-Hesse et al. 2003 ). 

The basic data for the galaxies and the complete log of 
the observations are listed in Table [T] and Table |2] respectively. 
5-band images of the galaxies are displayed in Figure T] IFS 
of III Zw 102 and Mrk 314 has been previously published by 
[ Garcia-Lorenzo et al.| ( |2Q08] ), who studied the central 33r6x29'.'4 
of III Zw 102 and W x 12'' of Mrk 314 using INTEGRAL at 
the William Herschel Telescope. In both cases, the much smaller 
field of view (FOV) of INTEGRAL allows only for a partial 
mapping of the starburst regions. 

Spectral data for the five galaxies listed in Ta ble \l\ were 
collected using the VIRUS-P IFU spectrograph ( [Hill et aT 



2008 ), working at the 2.7m Harlan J. Smith Telescope at the 
McDonald Observatory (TX, USA). VIRUS-P is the prototype 
instrument for VIRUS (= Visible Integral-field Replicable Unit 
Spectrograph), a massively replicated IFU made of 150 units like 
VIRUS-P, which will be installe d on the 9.2 m Hobby-Eberly 
Telescope ( [Hill et al.|20Q6l[2Q07l ). 

VIRUS-P (with the VP-1 IFU bundle used in these observa- 
tions) consists of a square array of 247 optical fibers (of which 
three were dead in the used configuration), which samples a 
109'' X 109'' FOV with a 0.3 filHng factor. Each fiber is circular 
and has a diameter of 4'.' 2 on the sky. We observed in the spectral 
range 3550-5850 A, with a spectral resolution (FWHM) of 5 
A (340 km s"^ at Bfi). 

Observations were carried out during two observing runs, in 
2007 October and in 2008 March. During the latter run, due to 
camera alignment issues, the spectra of five fibers fell off' the 
chip, and therefore only 239 fibers were usable. 
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Galaxy 


Other designations 


RA 


DEC 


Mb 


D 


Mb 


Morph 


Ay 






(2000) 


(2000) 




(Mpc) 


(mag) 




(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


II Zw 33 


UGCA 102, Mrk 1094 


05 10 48.1 


-02 40 54 


14.57^ 


37.4 


-18.29 


Ch 


0.288 


Haro 1 


NGC 2415, UGC 3930 


07 36 56.7 


35 14 31 


\2.6V 


52.9 


-21.00 


Ext 


0.117 


NGC 4670 


UGC 7930, Arp 163, Haro 9 


12 45 17.2 


27 07 32 


13.02^ 


22.4 


-18.73 


Nuc 


0.041 


Mrk314 


NGC 7468, UGC 12329 


23 02 59.2 


16 36 19 


13.78^ 


29.0 


-18.53 


Ch 


0.234 


IIIZw 102 


NGC 7625, UGC 12529, Arp 212 


23 20 30.1 


17 13 32 


12.50^ 


23.1 


-19.31 


Ext 


0.069 



Notes. Cols. (3) and (4): Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. 
Col. (5): ^-b and magnitudes: asymptotic photometry obtained by extrapolating the growth curves, and corrected for galactic extinction (Cairos] 
|et al.|200Ta| . The asymptotic magnitudes listed in Cairos et al. (2001a) were corrected for Galactic extinction fo llowing [B ur stein & Heiles| |7982| ); 
here they have been recomputed using the Schlegel et al. (1998) extinction values; integrated magnitudes from Gil de Paz et al.,(2003^ ), corrected 
for Galactic extinction. Col. (6): Distances, computed assuming a Hubble flow, with a Hubble constant Hq - 73 km s~^ Mpc , and taking into 
account the influence of the Virgo cluster, the Great Attractor and the Shapley supercluster, are taken from NED (http://nedwww.ipac.caltech.edu). 
Col. (7): Absolute magnitude in the B band, computed from the tabu lated B-band magnitu de and distance. Col. (8): morphological classification 
following [Cairos et"ar] { [2001^ . Col. (9): galactic extinction, from the |Schlafly & Finkbeiner|p011| recalibration, as listed in NED. 



We followed the same strategy in the two observing cam- 
paigns. Calibration frames consisted of bias frames, sky-flats 
and comparison lamps taken at the beginning and at the end of 
the night. Sky-flats are required to define and trace the fibers 
(apertures) on the detector as well as to perform the throughput 
correction. Spectra of emission-line lamps (Hg-Cd lamp in our 
case) are needed to carry out the wavelength calibration. Because 
VIRUS-P is an instrument with practically no flexures, compar- 
ison lamp spectra taken once during the night will suflftce, and 
no arc exposures interspersed with the science observations are 
required. 

As the filling factor of the IFU is 30%, three dithered expo- 
sures were taken in order to ensure a coverage of about 90% of 
the FOV. The standard observing sequence consists of a set of 
three exposures (of 20 to 30 minutes each) at each of the three 
positions (the telescope offsets are first 3V9 west, 7!\ south and 
then 3^75 east, 7! I south). 

Flux spectrophotometric standards were also observed, al- 
beit using a six position dither pattern, which guarantees the col- 
lection of (almost) their total flux. The spectrophotometric stan- 
dard stars BD-h25d4655, G 191-B2B and BD-h40d4032 were ob- 
served during the October run, and Feige 34 and BD-h33d2642 
in the March run. 

2.2. Data reduction 

The data have been processed using standard IRAlj^tasks. 

The first step in the data reduction was to remove the bias 
pedestal level, whose value was computed using the overscan 
regions in the science frames, and subtracted out. As the bias 
frames showed some smooth large-scale structure, we averaged 
all of them to obtain a master bias, to which we subtracted its 
mean value. This master bias was then subtracted from all the 
overscan-corrected frames. 

Next, apertures were defined and traced on the detector. The 
apertures need to be traced and defined in well exposed frames: 
sky-flats taken at the beginning or at the end of the night were 
used for that. We defined the apertures using the IRAF task 
apall\ this task first finds the centers of each fiber (the emission 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



peak) along the spatial axis at some specified position, and then 
the user is prompted for the size of the extraction window, which 
we set to 7 pixels, a width that, in our dataset, maxi mized the 
overall sign al-to-noise ratio (cross-talk is negligible, see| Murphy 
|et al.||20lT| ). The apertures were then traced by fitting a poly- 
nomial to the centroid along the dispersion axis. A fifth degree 
Legendre polynomial was found to provide good fits, with a typ- 
ical RMS of about 0.02 pixels. 

Once the apertures were defined and traced in the continuum 
frames, we used again apall to extract them in all the science 
frames. The extraction consists of summing, for each fiber, the 
pixels along the spatial direction into a one-dimensional spec- 
trum. The output is thus a "row-stacked spectrum" (RSS) of the 
244 spaxels (239 for the 2008 March run). 

We then performed the wavelength calibration. First, in the 
comparison spectra we identified several emission features of 
known wavelength in a reference fiber. Second, using identify, 
a polynomial was fitted along the dispersion direction; the stan- 
dard deviation (RMS) of the polynomial fit gives an estimate 
of the uncertainty in the wavelength calibration. We obtained a 
typical RMS of about 0.01 A by fitting a fifth degree Legendre 
polynomial. A total of 13 lines were fitted. Finally, with reiden- 
tify we identified the emission lines in all the remaining fibers of 
the arc frame, using the selected one as a reference. 

Subsequently we corrected for throughput (spaxel-to-spaxel 
variations) by using the sky flat exposures. This step was carried 
out with the task msrespl. 

In each spectrum, the sky background was computed by tak- 
ing the median of ~ 40 - 50 spaxels in the periphery of the 
FOV where the galaxy line and continuum emission is negli- 
gible, and subtracted out. The spectrum was then corrected for 
atmospheric transparency variations, gauged by means of the 
flux of the guide star in the guider saved frames, as explained 
in Blanc et al.|p00"9] ), and for atmospheric extinction (the guide 
star exposures were also used to estimate the seeing during the 
observations, see Table [2]). 

The three spectra for each dithering position were averaged 
together using the IRAF task imcombine, and then stacked to 
form a final spectrum of 732 rows (717 for the 2008 march ob- 
servations). 
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Table 2. Log of the Observations 



Galaxy 


Date 


Exp. time 


Airmass range 


Seeing 


Spatial scale 






(s) 




(arcsec) 


(pc arcsec" 


II Zw 33 


Oct 2007 


3600 


1.20-1.40 


2.2-2.4 


760 


Haro 1 


Mar 2008 


3600 


1.01-1.10 


2.5-3.2 


1080 


NGC 4670 


Mar 2008 


4800 


1.00-1.13 


2.0-3.0 


460 


Mrk314 


Oct 2007 


5400 


1.03-1.20 


1.8-2.0 


590 


III Zw 102 


Oct 2007 


4200 


1.03-1.25 


2.0-2.8 


470 




Fig. 1. 5-band images of the galaxies (in parenthesis the instrument/telescope at which the images were obtained): II Zw 33 (CAFOS 
at the 2.2 m CAHA,.Cair6s et al. 2QQlb) ; Haro 1 (Prime Focus at 3.5-m CAHA, ,Cair6s et al..2QQlbj; NGC 4670 (60 inch telescope on 
Mt. Palomar, |Girde Paz et al.|20Q3| ); Mrk 314 (prime Focus at the 3.5-m CAHA'^ ICaiFos et al.|2001bl ) and III Zw 102 (ALFOSC/NOT, 
unpublished). North is up, east to the left. The FOV is 109'' x 109'' (the FOV of the VIRUS-P instrument). All images are shown in 
logarithmic scale. 



2.3. Flux calibration 

Spectra of several spectrophotometric standard stars were ob- 
tained and reduced in the same way as the galaxy spectra. 

The only remarkable difference is that the standard stars 
were observed with a six-position dithering pattern, which en- 
sures that only about 3% of FOV is "lost" in inter-fibers gaps, 
while about two thirds of the area is covered by two overlapping 
fibers, as shown in [2] Combined with the PSF undersampling, 
this dithering pattern does not allow an accurate direct flux esti- 
mation by simple summation of all the relevant spaxels. 

If we summed together all the spaxels with sufficiently strong 
stellar continuum (typically about 30 spaxels in the row-stacked 
final spectrum), we would overestimate the actual stellar spec- 
trum, because of the dithering geometry. Moreover, the overes- 
timation factor is in principle wavelength-dependent, because of 
both the increase of the stellar PSF width towards bluer wave- 



lengths, and the shift of the light centroid due to diff'erential at- 
mospheric refraction. 

To better understand it, we can think of a hypothetical star 
with an extremely narrow PSF, whose blue centroid happens to 
fall on one of the inter-fibers gaps while the red centroid falls on 
a spatial location covered by two overlapping fibers: we would 
measure very little or zero blue flux, while the red flux would be 
overestimated by a factor of two. 

To assess the importance of this effect, we simulated the re- 
sponse of each fiber and did a fit to determine the PSF pa- 
rameters. The fiber response is calculated by a Monte Carlo in- 
tegration of the PSF over the circle area of the fiber. For the PSF, 
a Moffat profile is chosen to follow the instrument response to 
the outer wings, as shown in Fig. [3] Fit parameters were the total 
flux, source position and the PSF width a and exponent jS. The 
Moffat fit was done individually for each wavelength bin. 

We found a tight correlation between the Moffat function 
parameters and wavelength, as depicted in Fig. |4j also, the ob- 
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5x [arcsec] 



Fig. 2. Size and position of the fibers, centered on the position 
of the star (G 191-B2B). The fiber diameter is 4V2, the distance 
between adjacent fibers in the same exposure is 7! 5. The blue 
circle shows the computed FWHM of this star V!A). 
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Fig. 3. Measured and fitted intensity for the source G 191-B2B. 
The X axis shows the distance to the fitted source center, Y the 
intensity. Stars indicate the positions of the fiber centers. The 
curves are the best Moff'at (solid) and Gaussian (dotted) fits (both 
integrated over the fiber area). 



served shift of the PSF centroid agrees well with the shift com- 
puted theoretically from the airmass and the parallactic angle of 
the spectrophotometric standard stars spectra. 

However, the overestimation factor turns out to have only a 
very small (< 0.03 mag) dependence on wavelength (Fig. [5|, 
because of the relatively poor seeing, which spreads out the star 
flux over a large number of spaxels, and of having observed the 
standard stars at low airmasses. 

To derive the final flux calibration function, for each spec- 
trophotometric standard star we summed together all spaxels 
with sufficient strong signal in the stacked spectrum, and applied 
the overestimating factor correction. Then the IRAF tasks stan- 
dard and sensfunc were used to compute the calibration function, 
after combining together the data for the diff'erent spectrophoto- 
metric stars. 

By comparing the sensitivity curves for the diff'erent standard 
stars in the two observing runs, we can estimate that the relative 
uncertainty on the calibration factor is generally less than 2%, 
except at the blue end of the spectrum, where the uncertainty 
increases to about 7%. 

As the diameter of the fibers and the average seeing values 
are much larger than diff'erential atmospheric refraction (DAR) 
eff'ects, no correction for DAR was applied to the data. 

2.4. Emission-line fit 

2.4.1 . Recovering the pure emission-line spectrum 

The measurement of the emission-line properties (e.g. fluxes, 
equivalent widths, or velocity dispersion) provides a great deal 
of information in regard to the nebula physical conditions and 
the properties of the ionizing stellar population. However, and in 
particular when dealing with SF galaxies, the derivation of accu- 
rate emission-line fluxes and equivalent widths is not free from 



problems. At any given spatial position, the observed galaxy 
spectrum is the sum of the nebular and stellar emission; under- 
lying stellar absorption can under certain circumstances dramat- 
ically decrease the measured intensity and equivalent width of 
B aimer emission lines. 

Therefore, in order to obtain precise measurements of the 
emission-line parameters of the ionized gas, the flux emitted by 
the stellar population must be modeled and subtracted. In the 
case of IFU data, this process must be done for every spaxel, 
i.e., the technique chosen to model the stellar population has to 
be applied to each individual spectrum of the IFS mosaic. 

In order to model the stellar population, we us ed the penal- 
ized pixel-fitting (pPXF) methocj^ developed by Cappella ri &| 
Emsellenii] ( |2004| ); the pPXF software allows to simultaneously 
fit the optimal linear combination of stellar templates to the ob- 
served spectrum and to derive the stellar kinematics, using a 
maximum likelihood approach to suppress noise s olutions. We 
made use of the MILES stellar libraryp]( |Sanchez-Blazquez et al. 
|2006 ), which contains a total of 985 stars spanning a large range 
in stellar age and metallicity. 

As a first step, we produced, for every galaxy, a high signal- 
to-noise integrated spectrum, adding up all the spectra of the 
science spaxels. Using this spectrum we find the "optimal stel- 
lar template" for every object, i.e., the selected stars from the 
MILES library that produce the best fit. The fit to the stellar 
continuum was done after masking out the gas emission-lines 
expected in the spectral range under study, as well as residuals 
from skyline subtraction. The number of stars combined to cre- 
ate the "optimal stellar template" for each galaxy was as follows: 



^ URL: http://www-astro.physics.ox.ac.uk/ mxc/idl/ 
^ URL: http://miles.iac.es/ 
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Fig. 4. The Moffat a and /3 parameters are plotted as a function 
of the wavelength for three spectrophotometric standards. The 
bottom panels show the computed FWHM (the solid line repre- 
sents the FWHM oc relation) and the shift of the centroid 
with respect to a reference point at A = 5750 A. 



17 stars for II Zw 33, 20 stars for Haro 1, 22 stars for NGC 4670, 
20 stars for Mrk 314, and 25 stars for III Zwl02. 

After that, the pPXF method was run on the spectra of each 
individual spaxel, using the "optimal stellar template" and al- 
lowing the weights of each star to vary. The application of the 
"optimal stellar template" to each individual spectrum consider- 
ably speed up the fitting process. As the automatization of this 
process could lead to bad fittings, a careful individual inspection 
by eye was done. 

Once the stellar continuum has been successfully reproduced 
by a sum of templates for every spaxel, it was subtracted from 
the original spectrum to get the pure emission-line spectrum. 

We must remark here that the fit to the stellar continuum is 
not meant to provide a physical description of the stellar popu- 
lations making up the galaxy, a topic which is outside the scope 
of this paper, but that the purpose of this modeling is to obtain a 
good representation of the underlying stellar continuum in order 
to decouple it from the emission lines produced by the ionized 
gas, and to derive more accurate values of emission-line fluxes. 

The stellar continuum fit only works out in the central re- 
gions of the galaxies, while it fails, producing meaningless re- 
sults, in the outer regions where the signal-to-noise of the con- 
tinuum is too small. 
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Fig. 5. Ratio between the flux computed by a direct sum of the 
spectra of all spaxels with sufficiently good signal, and the flux 
derived by the Mofl'at fit to the data, as a function of wavelength, 
for three spectrophotometric standard stars. The direct integra- 
tion overestimates the stellar flux by a factor of about 5/3, which 
depends very little on wavelength. 



2.4.2. Fitting the lines 

Emission-line fluxes for every spaxel were measured on the pure 
emission-line spectrum. However, the line emission extends well 
beyond the galaxy region in which the pure emission-line spec- 
trum could be computed. 

To compute the fluxes in those spaxels outside such region, 
the original stacked spectra were used, fitting the continuum, 
typically 30-50 A on both sides, by a straight line. 

In order to measure the relevant parameters of the emis- 
sion lines (position, flux and width), they were fitted by sin- 
gle Gaussians. The fit procedure was carried out by using the 
the Levenberg-Marquardt non-linear fitting algorithm imple- 
mented by C. B. Markwardt in the mpfitexpr IDL library |^ 
The [Oiii] /1/14959, 5007 doublet was fitted imposing that both 
lines have the same redshift and width. Uncertainties on the fit 
output parameters were computed adopting a simplified poisso- 
nian model on the non-flux-calibrated spectra, where the readout 
noise was set equal to the rms of the nearby continuum, and the 
gain set to 1 . 

Criteria such as flux, error on flux, velocity and width were 
used to do a first automatic assessment of whether to accept or 
reject a fit. For instance, lines with too small (less than the in- 
strumental width) or too large widths were flagged as rejected, 
as well as lines with an error on the flux more than about 10% of 
the flux (the exact limits depending on the specific line and on 
the overall quality of the spectrum). Such criteria were comple- 
mented by a visual inspection of all fits, which led to override, 
in a few cases, the automated criteria decision (by accepting a fit 
flagged as rejected, or vice versa). 



2.5. Creating the 2D maps 

The emission-line fit procedure gives, for each line and for 
each line parameter (for instance flux), a table with the fiber 
ID number, its coordinates, the measured value and the accep- 
tance/rejection flag. 

The gaussian fit parameters for the continuum- subtracted 
and the original spectra were merged together. For the Hy and 



URL: http://purl.com/net/mpfit 
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the US B aimer lines, due to their strong absorption components, 
only the results from the continuum- subtracted spectra were 
used. 

Line ratios were computed by simply dividing the fluxes of 
the corresponding lines. 

To display all these maps, we explored several options. One 
is to interpolate the data on a r egularly s paced grid, by means of 
an algorithm such as the Renka & Cline| ( 1984 ) two-dimensional 
interpolation method. However, our own experience is that the 
interpolation can produce spurious oscillations, especially over 
non-contiguous spaxels, and gives the false impression of a spa- 
tial resolution higher than the actual value. 

An alternative is to represent each spaxel with measured data 
with a circle with the fiber's diameter at its spatial position. 
However, the irregular inter-fiber gaps (due to the fact that the 
actual position of the fibers is slightly and randomly off'set from 
the nominal one) distract the eye of the viewer and interfere with 
a comfortable and eff'ective visual assessment of the global char- 
acteristics of the map. 

The solution we adopted was to depict each fiber by means 
of an hexagon centered on the nominal (not the actual) position 
of each fiber (the mean diff'erence between nominal and actual 
positions is about Off 4 with a maximum of IV 8). This we found 
to be the best compromise between a truthful representation of 
the data (as no rebinning or interpolation was done: the only ma- 
nipulation was to slightly change the fibers positions and shape) 
and a map which is both pleasant to the eye and easy to appraise 
and interpret. 



3. Results 

3.1. Emission-line and continuum intensity maps 

IFS provides a simultaneous mapping of the light emitted by a 
galaxy in a wide spectral range. Accordingly, from one single 
exposure we can retrieve maps at a given wavelength or maps 
equivalent to narrow- or broad-band images within the whole 
observed range. 

Following the steps described 



m Section 2.5 we built 
emission-line flux maps for the brightest observed emission- 
lines, namely [O ii] A3121, Hy, UJ3, and [O iii] A5001; these maps 
are displayed in Figure |6] (rows 2 to 5). 

Continuum maps were obtained by summing the flux within 
specific wavelength intervals, selected so as to avoid emission 
lines or strong residuals from the sky spectrum subtraction. 
Continuum maps in the spectral interval 4500-4700 A are shown 
in Figure |6] (row 1): this is the wavelength range where the con- 
tinuum signal is stronger and better portrays the galaxy outer 
regions (Wolf-Rayet features are not an issue, as they are only 
visible in just a few central spaxels in two objects, see Sect.|4]). 

The emission-line maps depict the interstellar gas ionized by 
massive stars, while the continuum maps mainly trace the stel- 
lar light from the galaxy. The five galaxies display an irregular 
morphology in emission-lines; in II Zw 33 and Mrk 314, the SF 
regions are aligned to form a chain, while in Haro 1 , NGC 4670 
and III Zw 102 the ionized gas departs from the inner regions, 
giving rise to a complex pattern, where the presence of several 
bubbles and filaments is remarkable. 

On the other hand, the continuum morphology is regular in 
the five objects; (the jagged appearance of the Haro 1 continuum 
map is due to the presence of a few foreground stars, which have 
been masked out, and a foreground galaxy at the west). In all the 
galaxies, the continuum peak is located roughly at the center of 



the outer isophotes, and its position does not depend on which 
specific wavelength range is used to build the continuum maps. 

3.2. Line ratio maps 

The first row of Figure |7] displays the [O iii] A5001/UJ3 maps 
for the five studied galaxies. This ratio essentially probes the 
hardness of the stellar radiation field within the nebula, and so 
high values of [O iii] A5001/UJ3 are predicted when the ioniza- 
tion is produced predominantly by UV photons from OB stars 
( [McCall et al. 19 85 ). Howe ver it is also sensitive to other fac- 
tors as e.g. the metallicity (Searle||1971| ), t he distance from the 



ionizing source to the [ O m] emitting layer (|Stasinska|2009 



or 

the presence of shocks ( [Dopita & Sutherland! 1995| ). Therefore 
the information provided by [O iii] A5001/UJ3 is limited unless 
additional data regarding these other parameters are added. 

For II Zw 33 and Mrk 314, the excitation maps roughly 
mimic the emission-line maps, the highest [O iii] A5001/UJ3 val- 
ues (> 3), spatially coinciding with the highest peak in emission- 
lines, as usual in systems photoionized by hot stars. 

Both Haro 1 and III Zw 102 show more complex 
[O III] A5001/UJ3 patterns, with the highest values of the excita- 
tion localized in the galaxy periphery. This behavior, frequently 
found in spiral galaxies, has been interpreted as an abundance 
gradient ( Smith 1975). Interestingly, Haro 1 and III Zw 102 are 
the two most luminous objects in the sample and both relatively 
far from the dwarf galaxy regime. 

NGC 4670 shows a somewhat complex excitation morphol- 
ogy, with several [Om] ASOOl/HjS peaks, located in the central 
giant H ii region and at the outer edge of the two larger filaments. 
For all five galaxies, the excitation levels are consistent with H ii- 
like ionization. 

Interstellar extinction can be probed by comparing the ob- 
served ratios of hydrogen recombination lines with their theoret- 
ical values ( [Osterbrock & Ferland|2006| ). In the optical domain, 
extinction maps are usually derived from the ratio Ha/Hfi as it 
involves the strongest lines, and therefore, the easiest to measure. 
In our case, lacking Ha, we have derived the extinction from the 
HyS/Hy line ratio; the extent of the extincti on m aps is dictated by 
the extent of the Hy maps (as said in Sect. 2^ we measured Hy 
only in the continuum- subtracted spectra. 

The extinction maps are shown in the second row of Fig. |7] 
using the Ay parameter. We adopted the Case B, low density 
limit, T = 10000 K approximation, where the theoretical H/5/Hy 
ratio is 2.15 (Osterbrock & Ferland 2006). The Cardelli extinc- 
tion law was used to derive Ay from the observed UjS/Hy value. 

3.3. Ionized gas and stellar kinematics 

We study the kinematics of the ionized gas using the 
[Om] A5001 emission-line, except for III Zw 102 for which we 
use UJ3, much stronger than [O iii]. Due to the low spectral reso- 
lution of our data (about 5 A FWHM), no reliable velocity dis- 
persion measurements could be obtained. 

The velocity fields are shown in Figure[7](row three); in these 
maps red and blue colors represent, respectively, redshifted and 
blueshifted regions with respect to the recessional velocity of 
each BCG. To estimate the errors on the velocity data, we made 
a combined plot (all galaxies together) of the uncertainty on the 
[O III] A5001 emission-line centroid versus the logarithm of its 
flux. For log(flux) > 4 the error is < 5 km s"\ in agreement with 
the mean uncertainty estimated from the scatter in the velocity 
maps built on the 5577 A skyline. The error increases to about 
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Fig. 6. Observed intensity maps (not reddening-corrected). First row: continuum emission between 4500 and 4700 A (a spectral 
region free from gas emission). Rows 2 to 5: emission-line flux maps derived from Gaussian fits to each line and spaxel (see the text 
for details). Axis units are arcseconds; north is up, east to the left. All images are shown in logarithmic scale. Flux units are 10"^^ 
ergs cm"^ s"^ Spaxels where the blue Wolf-Rayet bump is visible are marked by stars in the [O iii] map. 



10-15 km s-i at log(flux) = 3, to 25-30 km s'^ at log(flux) = 2, 
and reaches or exceeds 40 km s"^ at fainter flux levels. 

The pPXF program provides, besides the best fitting stellar 
template, the radial velocity and associated uncertainty for each 
fitted spaxel, allowing us to produce the galaxy's stellar velocity 
field (given in Figure |7j row four). 

The uncertainty on the stellar radial velocity varies between 
^10 km s"^ in the nuclear regions and ^ 40 km s"^ in the outer 
spaxels. 



3.4. Integrated spectroscopy 

For every galaxy we have also produced its integrated and nu- 
clear spectrum. In order to generate the integrated spectrum we 
included only those spaxels in which the stellar contribution was 
successfully modeled as a sum of templates (see Section [24| ); 
in this way we ensure a better correction for the contribution 
of the underlying stellar population and, therefore, a more ac- 
curate determination of the emission-lines fluxes. The integrated 
spectrum was produced by summing together the individual pure 
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Fig. 7. First row: Observed (not reddening-corrected) ionization ratio [O iii] A5001/11/3 for the five galaxies. Second row: Interstellar 
extinction Ay map, computed from the HjS/Hy ratio. Row three: Velocity field (in km s"^) of the ionized gas [O iii] /15007 emission 
line, except for III Zw 102 for which the H/S velocity field is shown. Row four: Stellar velocity field. Axis units are arcseconds. The 
solid line in the gas velocity map of NGC 4670 is the isovelocity contour at V = 1050 km s"^ (see text for details). 



emission-line spectra of the selected spaxels. A corresponding 
continuum spectrum was obtained by summing the best-fit con- 
tinuum spectra for the same spaxels. 

On the same line, we also obtained line-emission and best-fit 
continuum spectra for the central region by summing over 6-8 
spaxels around the continuum peak (for short, we shall refer to 
them as "nuclear spectra"). 

Figure [5] displays the integrated spectra for the sample ob- 
jects. All objects but III Zw 102 display a typical BCG spec- 
trum, in which strong optical emission-lines are superposed on 
an almost flat continuum. On the contrary. III Zw 102's spec- 
tral shape, with a high and red continuum, resembles more those 
typical of an early type galaxy, rather than those of BCGs. 
All objects display prominent absorption features, including the 
higher-order Balmer lines, H^, Hy, HJS, Can H and K lines, and 
the G-band at 4304 A. 

For both the integrated and nuclear spectra, the flux of the 
visible emission-lines was measured on the pure emission-line 
spectrum using the Gaussian profile fitting option in the IRAF 



task splot. In order to compute equivalent widths, the contin- 
uum level was measured on the best-fit continuum spectrum. For 
Balmer lines, the continuum level was determined by interpolat- 
ing over the absorption feature using two continuum windows 
on either side. The interstellar extinction coefficient, C(tlj3), was 
computed through a least- square fit to the Balmer decrement, 
using the UJ3, US and Hy lines. 

Reddening-corrected intensity ratios and equivalent widths 
of the integrated and nuclear spectra are listed in Table [3] The 
listed uncertainties on the line ratios and the UJ3 fluxes are de- 
rived by summing in quadrature measurement errors, uncertain- 
ties on the flux calibration (Sect. 2.3 ) and the uncertainty on the 
computed C(UJ3). 



3.4.1 . Heavy element abundances 

In the two galaxies Mrk 314 and II Zw 33, where the 
[O III] /14363 auroral line could be measured, electron temper- 
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Table 3. Reddening corrected line ratios, normalized to UJ3, for the sample of galaxies 



A 


Ion 




II Zw 33 




Haro 1 






Integrated 


Nucleus 


Integrated 


Nucleus 






Fa 


-Wa 


Fa 


-Wa 


Fa -Wa 


Fa -Wa 


3727 


[On] 


3.312 ±0.284 


58.94 ± 1.43 


3.125 ±0.254 


48.26 ± 0.70 


2.419 ±0.217 24.25 ±0.32 


2.106 ±0.190 19.55 ±0.15 


3835 


H9 


0.055 ± 0.006 


0.80 ± 0.05 


0.041 ± 0.004 


0.50 ± 0.04 








3868 


[Nelll] 


0.259 ±0.018 


3.53 ± 0.06 


0.236 ± 0.015 


2.66 ± 0.04 


0.062 ± 0.006 0.44 ± 0.02 


0.040 ± 0.004 0.27 ± 0.01 


3889 


Hel+H8 


0.166 ±0.012 


2.39 ± 0.04 


0.168 ±0.011 


1.95 ±0.03 


0.132 ±0.009 0.94 ±0.02 


0.126 ±0.009 0.85 ±0.01 


3968 


[Nelll] 


0.228 ± 0.014 


3.21 ± 0.04 


0.209 ± 0.012 


2.47 ± 0.04 


— — 


— — 


4101 




0.256 ± 0.012 


3.68 ± 0.03 


0.256 ±0.012 


3.02 ± 0.06 


0.244 ±0.013 1.84 ±0.02 


0.239 ±0.012 1.73 ±0.01 


4340 


Hr 


0.479 ± 0.017 


7.82 ± 0.07 


0.476 ± 0.015 


6.26 ± 0.06 


0.469 ± 0.018 3.96 ± 0.01 


0.469 ±0.018 3.88 ±0.01 


4363 


[OUT] 


0.032 ± 0.003 


0.51 ±0.03 


0.024 ± 0.002 


0.31 ±0.03 






4861 




1.000 


19.57 ±0.13 
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15.57 ±0.09 


1.000 10.49 ±0.05 
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4959 


[OIII] 


0.903 ± 0.023 


18.66 ±0.10 
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2.554 ± 0.064 
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0.532 ± 0.013 6.22 ± 0.02 






/(I^) = 1061.6 ±18.5 


I(B/3) = 489.6 ± 8.5 


I(BJ3) = 4185.9 ± 72.8 


I(UJ3) = 1532.0 ± 26.5 






F(U/3) = 2090.3 ± 283.3 


F(U/3) = 760.0 ± 80.4 


F(U/3) = 9800.7 ± 1628.1 


F(U/3) = 3600.5 ± 597.4 






C(HyS) = 0.294 ± 0.058 


C(U/3) = 0.191 ±0.045 


C(H/3) = 0.369 ± 0.072 


C(U/3) = 0.371 ±0.072 






Ay - 0.630 ±0.125 


Ay = 0.409 ± 0.097 


Ay = 0.791 ±0.154 


Ay = 0.794 ±0.154 
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= 24.3 kpc^ 
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' = 4.5 kpc2 
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[On] 


3.051 ±0.280 
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2.487 ± 0.200 


44.16 ±0.30 


3.255 ± 0.270 50.66 ± 0.92 


2.859 ± 0.264 54.09 ± 0.60 


3835 


H9 


0.048 ± 0.004 


0.89 ± 0.03 


0.044 ± 0.004 


0.71 ±0.01 


0.053 ± 0.005 0.97 ± 0.07 


0.037 ± 0.003 0.78 ± 0.02 


3868 
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3.32 ± 0.03 
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3.37 ±0.01 


0.296 ± 0.019 3.58 ± 0.06 


0.297 ± 0.022 4.24 ± 0.02 


3889 


Hel+H8 


0.187 ±0.014 


3.18 ±0.02 


0.182 ±0.011 


2.84 ± 0.01 


0.165 ±0.010 2.91 ±0.08 


0.148 ±0.010 2.12 ±0.02 
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[Nelll] 


0.217 ±0.013 
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0.215 ±0.010 


3.55 ±0.01 


0.245 ±0.013 3.34 ±0.03 
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0.469 ± 0.017 7.88 ± 0.04 
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1.000 


26.03 ±0.14 
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ro Tin 


2.223 ± 0.056 


54.01 ± 0.23 


2.500 ±0.061 


72.26 ± 0.23 


3.053 ± 0.075 59.08 ± 0.55 


3.354 ± 0.085 77.07 ± 0.63 






7(1^) = 6636.0 ± 116.1 


I(Ufi) = 3101.8 ±54.3 


I(U/3) = 1491.8 ±26.2 


I(U/3) = 992.0 ± 17.3 






F(Ufi) = 10489.3 ± 1598.2 


F(U/3) = 5565.2 ± 382.6 


F(iy3) = 3614.7 ± 394.0 


F(Hfi) = 1976.5 ± 330.3 






C(I^) = 0.199 ±0.066 


C(Hfi) = 0.254 ± 0.029 


C(Ufi) = 0.384 ± 0.047 


C(Ufi) = 0.299 ± 0.072 






Ay = 0.426 ±0.141 


Ay = 0.543 ± 0.062 


Ay = 0.823 ±0.100 


Ay = 0.641 ±0.155 
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0.237 ± 0.015 


1.01 ±0.01 


0.233 ±0.012 
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0.469 ± 0.019 


2.20 ± 0.02 


0.469 ±0.018 
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1.000 


4.93 ± 0.02 


1.000 
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[OIII] 


0.112 ±0.004 


0.58 ± 0.01 


0.086 ± 0.003 


0.76 ±0.01 
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[OIII] 


0.319 ±0.008 


1.77 ±0.01 


0.230 ± 0.006 


2.18 ±0.01 










I(B/3) = 1652.9 ± 29.2 


I(B/3) = 809.2 ± 14.2 










F(I^) = 6031.4 ± 1075.0 


F(BJ3) = 3647.1 ±628.9 










C(HyS) = 0.562 ± 0.077 


C(Ufi) = 0.654 ± 0.075 
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Ay = 1.399 ±0.160 
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Notes. Reddening-corrected line fluxes, normalized to F(H/?) = 1. The observed Hfi flux, I(tlj3), the reddening-corrected H/? flux, F{llj3) (both 
xlO~^^ ergs cm~^ s~^), the reddening coefficient, C(H/?), and Ay (derived as 2.14 x C(H/?)) are listed for each region. Also shown are the area 
covered by the integrated spectra, and the corresponding (luminosity weighted) average stellar and gas velocities. The quoted uncertainties account 
for measurement, flux-calibration and reddening coefficient errors. 



atures and oxygen and neon abundances were derived using the 
"direct" method. 



First, we computed the electron temperature re(OIII) from 
the [Oiii] /l 4363/[Qiii] /14 959 + /15007 ratio using the equations 
provided by AUer ( |1984 ). As in our spectra there are no useful 



lines to obtain the electron density, A^e, we assumed A^e 
cm"^. 



100 



The metal abundances were derived using the revised expres- 
sions in Iz otov et al. ( 2006b), which take into account the latest 
atomic data and incorporate an appropriate grid of photoioniza- 
tion models with state-of-the-art model atmospheres. We used 
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Fig. 8. Integrated spectra of the sample galaxies, obtained by 
summing together the spectra of all those spaxels in whic h th e 
stellar continuum could be fitted by a template (see § |2.4| ). 
Spectra are plotted in logarithmic flux scale, and have been 
shifted vertically. Major tick-marks on the y axis are spaced 1 
dex apart. Nuclear spectra do not exhibit significantly distinct 
features from the integrated spectra and are not shown. 



their equations (3) and (5) to compute the 0^/H, and 
Ne^^ abundances, after computing 7^(011) using equation (14). 
The total oxygen abundance was obtained by simply summing 
the abundances of and O^^, as the fraction of the O^^^ ion 
is generally negligible. The correction factor for Ne, /CF(Ne^^), 
was computed using equation (19). (Using the Ne^^/0^^ ratio 
instead of the ICF we would get a Ne abundance about 0.2 dex 
larger.) 

For the remaining three galaxies, where [O iii] /14363 could 
not be measured, we had to rely on the so-called "strong-line 
methods", in order to derive the oxygen abundance. We adopted 
the revised calibratio n of the P -method ( |Pilyugin|[2000| |2001| ), 
provided in Pilyugin & Thuan| (2QQ5 ). 

Electron temperatures and abundances are listed in Table |4] 
For II Zw 33 and Mrk 314, we provide only the abundances de- 
rived from the direct method; since both galaxies lie in the tran- 
sition zone 8.0 < 12 -h log(0/H) < 8.3, the Pilyugin & Thuan 
( |2005| ) calibration does not give reliable abundance values. 



3.4.2. Comparison with literature data 

To assess the quality of our measurements, we compared our 
integrated fluxes and equivalent widths with the integrated data 
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Fig. 9. Upper panel: ratio between the fluxes computed by di- 
rect summation of our integrated spectra, after correction for 
foreground galactic extinction, and the (g alactic extinction cor- 
rected) integrated fluxes published by Moust akas & Kennicutt 
(2006). Lower panel: ratio between equivalent width measure- 
ments. To avoid overlapping of error bars (which include the 
uncertainties on both datasets), datapoints have been slightly 
shifted from their nominal position on the X axis. 



published by [Moustakas & Kennicutt|p006| ). Four of our objects 
are in common with their sample, the exception being II Zw 33. 

For this comparison, our measurements have been corrected 
for foreground galactic extinction only. The results of the com- 
parison are displayed in Figure [9l Taking into account the difl'er- 
ent integration area used in the [Moustakas & Kennicutt| p006 
spectra (a 30'' x 50'' aperture) and our own spectra, the compar- 
ison looks fair. 



4. Discussion 

4.1. IIZw33(=Mrk 1094) 



II Zw 33 is a well- studied object, cataloged as BCG in [Thuan & 
Martin (1981). Optical imag ing and surface photometry has been 
presented in several works ([Loose & Thuan|1986 Cairos et al. 
[2001a|bl [Gil de Paz et al. [[2003 [ [Gil de Paz & Madore[[2005[ ). 



The presence of a central SF bar and the striking twistin g of 
the isophote s have been interpreted as interaction signs (Mendez| 
[eFaLlP^99l ). II Zw 33 is also a well kno wn Wolf-Ray et (WR) 
gal axy, included in the WR galaxy lists by [Vacca & Conti[ ( [l992 1 
and [Schaerer et al.[ ( [T999l ). 

Maps of II Zw 33 are displayed in Figures |6]and|7] At the 
II Zw 33 distance, the 4'.' 2 fiber diameter provides a spatial sam- 
pling of 760 pc per spaxel. 

The line emission is concentrated along a central bar-like 
structure, with no clear peak, while the stellar emission extends 
further out from the central SF regions, evidencing the presence 
of a much more extended, more evolved stellar population. The 
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Table 4. Electron temperatures and heavy metal abundances. 



Value 


II Zw 33 
Integrated Nuclear 


Haro 1 
Integrated Nuclear 


NGC 4670 
Integrated Nuclear 


Mrk314 
Integrated Nuclear 


III Zw 102 
Integrated Nuclear 


re(oiii) (K) 

TeCOII) (K) 
12 + log (0/H)i 
12 + log (0/H)2 
log(Ne/0) 


12290 + 330 11270 + 360 
12080 + 290 11160 + 350 
8.06 + 0.03 8.16 + 0.04 

-0.85 + 0.05 -0.85 + 0.06 


8.36 8.40 


8.29 8.37 


12330 + 250 12090 + 210 
12110 + 210 11910+190 
8.08 + 0.03 8.10 + 0.03 

-0.83 + 0.04 -0.82 + 0.04 


8.44 8.49 



Notes. Oxygen abundances computed using the method. Oxygen abundances derived from Pilyugin & Thuan||2005| ); uncertainties are 
about 0.1 dex. 



continuum map displays a well defined peak, located roughly at 
the center of the galaxy. 

The excitation map ([Oni]/HyS) traces the regions of star- 
formation, as expected for objects photoionized by stars. The 
peak in [O in]/Hy g is located at the southern tip of the bar (knot 
'A' as labeled in |Mendez et al.|19"9 9); their value of 5.5 agrees 
well with the peak value we find, ^4.5 (the small discrepancy is 
likely due to the diff'erence in spatial resolution). 

The extinction is relatively uniform in the central region, 
with Ay values around 0.4; these values are si gnificantly smaller 
than those published by |Mendez et al. (1999)). 

The [O in] velocity map does not show a clear rotation pat- 
tern, although the northwest side seems to have a slightly higher 
radial velocity as a whole. The ionized gas kinematics is in good 
qualitative agreem ent with the H I kinematical map published by 
( [Walter et al.|199 7, see their Fig. 6). The stellar velocity field is 
not as extended as required to make a sensible comparison, and 
is also much noisier. 

We derived an oxygen abundance of 12 -h log (0/H) = 8.06, 
in good agreement with the value of 12 -h log (0/H) = 7.99 re- 
ported by Shi et al.|p005| ), who also used the method. 

The Wolf-Rayet bump around 4600-4680 A is visible in two 
spaxels, marked by a star in the [O iii] flux maps in Figure |6] 

4.2. Haro 1 (=NGC2415) 



Haro 1 is classified as BCG by Gordon & Gottesman ( |1981 |) and 
Klein et aL](|1984|), and is included in the H ii galaxies sample 
studied by 'Deeg et al.| ( |1997 ). It has been defi ned as a paired 
galaxy with little (if any) tidal disturbance ( Sch neider & Salpeter] 
|1992| ). Its companion is UGC 3937, with a redshift diff'erence of 
only 210 km s"^ and a projected distance of 350 kpc (23 arcmin 
at Haro I's dista nce). Broad-band surfac e brightness photometry 
was presented in [Cairos et al.| (|2001a'b'). 

Maps are shown in Figures [6]and[7](column 2). At the Haro 1 
distance, the spatial resolution is about 1.1 kpc per spaxel. 

Haro 1 displays outer regular isophotes in the continuum 
maps. (The jagged appearance is the result of masking several 
foreground stars; for the sake of clarity, a smaller star, located at 
^ 30'' east from the center, has been left, as well as the back- 
ground galaxy on the west-side). The continuum emission peak 
is placed roughly at the center of the isophotes. The emission- 
line maps display a more irregular pattern: a considerably large 
filament extends out to the north-west, and the emission peak 
seems to be slightly off'set to the west with respect to the contin- 
uum peak. 

The [Oiii]/HyS ratio map has a shallow minimum, 0.40, in 
the galaxy's central region (where the emission peaks), and in- 
creases slightly outwards, reaching its maximum at the position 
of the extension detected in emission-line maps. Radial changes 



in the [Oin]/HyS ratio have been found in spiral galaxies, and 
have been inte rpreted as variation in the heavy element abun- 
dances ( Searle|[T97TJ [Smith[p775] ). Interestingly, Haro 1 is the 
most luminous object in our samp le, and it has been a lso classi- 
fied as a Sm galaxy (HyperLedap] [Paturel et al.|2003 ). 



The extinction map shows little or no spatial variations, with 
an average extinction Ay ~ 0.90. 

Haro I's kinematics show an overall regular rotation about 
an axis oriented southeast-northwest, with pretty much the same 
pattern and amplitude (about 160 km s"^) both in the gas and in 
the stellar velocity maps. 

The value derived for the metallicity, 12 -h log (0/H) = 8.36, 
is slightly lo wer that the value of 12 -h log (0/H) = 8.52 derived 
in Shi et al.[(2005j , which used the P-method (as calibrated in 
yugin|2001| ). 



4.3. NGC 4670 (=Haro 9, Arp 163) 

NG C 4670 has been classified as a blu e amorphous galaxy 
(Hu nter et al.[[1994[ Marlowe et al.[[1997| ), as a peculiar, barred 
so/a galaxy ([de Va ucouleu rs et al.|199T] ) and as a BCG (Gil de| 



Paz et al.][20b3| ). It is listed as a peculiar galaxy with "diff'use 
c ounter tails" in Arp's c atalogue ( Arp 1966[), and as aWR galaxy 
in [Mas-Hesse & Kunth[ ( |1999| ) and |Schaerer et al.[ ( [T999l ). Broad- 
band photometry and Ha imaging have been pro vided in [Gil de 
[Paz et al.[ ( [2Q03] ), [Gil de Paz & Madore[ ( [2005] ) and |Marlowe etaL 
(|1997|). It is also cataloged as a WR galaxy . 

Maps are shown in Figures [6] and [7] (column 3). At the 
NGC 4670 distance, the VIRUS-P IFU 4^2 fiber diameter corre- 
sponds to about 460 pc. 

The continuum map shows an outer regular morphology, 
with elliptical isophotes; the bulge, extremely bright, is clearly 
elongated on an east- west axis. The continuum peak appears to 
be a little displaced to the east, in agreement with the 5-band 
image (Figure [T]). In emission-lines, the galaxy displays a rather 
complex pattern, with a central, roughly circular and very bright 
Hn region, surrounded by an extended filamentary envelope, 
which has several bubbles, loops and extensions. Remarkable 
are the two biggest filaments, which expand to the southwest and 
northeast, extending up to about 7.5 and 6.5 kpc respectively. 

The ionization map displays an intriguing structure. The 
[O in]/HyS ratio clearly shows higher values in the central region, 
with a peak value ~ 3, and at the outer edge of the two larger ex- 
tensions. Highest values of [Oin]/HyS in the central knot, where 
the line-emission intensity also peaks, are consistent with pho- 
toionization produced by hot stars; it is not clear however what 
the ionizing mechanism acting on the filaments is. Although 
the morphology suggests a shock-like excitation mechanism, no 
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conclusions can be drawn without the measurement of additional 
diagnostic lines. The excitation values we found agree very well 
with those reported in Hunte r et al.| (p^994), [Oiii]/HyS = 2.97, 
derived for the spectra positioned on the center of the galaxy. 

The extinction map shows a roughly flat behavior in the inner 
region, with Ay around 0.40, and seems to be dominated by noise 
in the outer parts. 

Hunter et al. ( |1994 )), analyzing VLA data, report an "ordered 



NGC 4670 



rotation east- west, and the motion appears to be that of a disk at 

the same position angle as the optical galaxy ( 100 deg)". Our 

velocity maps show a more complex picture. At first sight a rota- 
tion pattern, though somewhat irregular, is visible in the gaseous 
component, with an amplitude of about 150 km s"\ around an 
axis oriented roughly northwest-southeast; the velocity gradient 
seems to be a little steeper in the northeast side of the galaxy 
compared to the southwest part. The stellar velocity field is over- 
all flat, and broadly agrees with the gas kinematics in the same 
area, except for the slight velocity rise toward its east edge, visi- 
ble in the [O iii] radial velocity map but not in the stellar one. 

The behavior of the gas and stellar kinematics can be better 
perceived in Fig 10 where we plot a velocity profile along a line 
with PA = 50 deg, where the gas distribution extends the most 
from the galaxy cent er, and a velocity profi le at PA = -100 deg, 
the value adopted byl Hunter et al.| ( p"994| ). Both graphs show a 
somewhat irregular rise of the gas velocity in a central region 
of about 40'' in size, and are flat, though noisier, outside it. The 
stellar velocity profile shows that the stars co-rotate with the gas, 
but with a lower velocity gradient, and seems to flatten out at 
somewhat smaller radial distance than the gas. Furthermore, the 
isovelocity contour at V = 1050 km s"^ in the [Oiii] velocity 
field is clearly S-shaped, a characteristic typical of barred galax- 
ies ( [Peterson et al.|1978] ). 

Thus, the kinematics of NGC 4670 is unlike the one of a 
rotating disk, and resembles more the kinematics of a barred 
galaxy, in line with its morphological classification, SB(s)0/aP?, 
listed in |de Vaucouleurs et al.| ( |199r] ). 

The derived metallicity, 12 -\- log (0/H) = 8.29, is sHghtly 
lower but consistent with the value of 12 -h log (0/H) = 8.4 



reported in ^Mas-Hesse & Kunth| ( |1999| ). 

The Wolf-Rayet bump around 4600-4680 A is visible in 
three spaxels, marked by a star in the [Oiii] flux maps in 
Figure [6] 

4.4. Mrk 314 (=NGC 7468) 

Mrk 314 has been cataloged as a polar-ring galaxy candidate 
dWhitmore et al.||1990| [van Driel et al.|[200n| ) and is included 
in the I Mazzarella & Boroson ( 1993 ) catalogue of multiple nu- 
clei galaxies. Su rface photometry in the optical and in the NIR 
was presente d in [Cairos et al. ( 2001a|b[ ), Caon et aL] ( |2005j ) and 
[Noeske et al.1 ( [20Q5| ). 

The maps of Mrk 314 are shown in Figures [6] and [7] (col- 
umn 4). The spatial sampling is about 590 pc per spaxel. 

The galaxy displays a similar morphology in both emission- 
line and continuum maps, with the SF knots aligned along the 
northeast southwest direction. In the continuum map, a more 
extended stellar component, with elliptical isophotes, surrounds 
the SF region. The brighter knot is located roughly at the cen- 
ter of the outer isophotes, and a tail-like feature departs in the 
southwest direction; a second filament extends out to the north- 
west. In the emission-line maps, the emission is concentrated in 
a central SF bar, the emission peak being slightly displaced from 
the continuum peak. 
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Fig. 10. Upper panel: Kinematical profile of NGC 4670 along a 
spatial cut through the reference point (0.0,8.0) and PA = 50 deg. 
The X axis is the distance of the projection of each spaxel on the 
above line to the reference point; only spaxels whose distance 
to the line is less than 6'' have been selected and plotted. Lower 
panel: Kinematical profile through the same point as above and 
PA = -100 deg. Open points are [Oiii] velocities, solid points 
stellar velocities. 



The peak of the [O iii]/HyS ratio is spatially co-located with 
the line emission maximum, as expected for objects ionized by 
OB hot stars; the excitation peak value of ^ 4.0 agrees well with 
the value found by jGarcia-Lorenzo et al.| ( [2008| ). 

The extinction map does not show any definite pattern, but 
a more or less homogeneous distribution with somewhat higher 
values at the south, and a Ay value around 0.7-0.8 in the central 
part; this result is consisten t with the Ha/ll/S map presented in 
[Garcia-Lorenzo et aLl ( [2008| ), measured in a smaller FOV (16'' x 
12''). 

The ionized gas velocity field displays a complicated pat- 
tern. Around the main body of the galaxy, the velocity field is 
similar to a rotating disk, with the rotation axis oriented east- 
west, the northern side approaching and the southern one reced- 
ing. This rotation is however perturbed in the central galaxy re- 
gions (within 10"), where the velocity field shows an S- shape 
distortion, indicative of the presence of a second rotating system. 
South of the main body (this area coincides with the tail-like fea- 
ture), the ionized gas is not rotating with the main-body galaxy 
disk. 

An extensively kinematic analysis of this galaxy has been 
done by j Shalyapina et al.[ p004[ ), who found that Mrk 314 
harbors two kinematics subsystems of ionized gas: the main 
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gaseous disk, and a smaller inner disk oriented orthogonally to 
it. Our coarse spatial resolution does not allow us to either con- 
firm or disprove the existence this inner disk, but the kinematics 
of the main disk agree well with the data shown by Shalyapina 
et al. ( 2Q04| ): see for instance the peak in the velocity profile at 
^ 30'' south of the nucleus, and then the sharp decrease (their 
Fig. 2a). 

The stellar velocity field looks more irregular, though there 
seems to be some rotation along an axis at about PA =~ 45 deg, 
with a total amplitude of about 70-80 km s"^ 

The oxygen abundance we found in the integrated spectrum, 
12 + log (0/H) = 8 .08, is somewhat lower than the one reported 
by |Shi et al.| ( [2QQ5] ), 12 + log (0/H) = 8. 26 (derived using also 
the Tq method), and the ones reported by [Garcia-Lorenzo et al.| 
( 2008| ), which range spatially from 8.14 to 8.35 (using the em- 
pirical calibrator N2, |Denicol6 et al.|2002| ). 



4.5. \\\Zw 102 (=NGC 7625, Arp 212) 

III Zw 102 is a galaxy well studied in a large range of 
wavelengths, but still a puzzling system. It is included in 
the Vorontsov-Velyaminov catalogue of interacting galaxies 
( [Vorontsov-V elyaminov 1959) and in Arp's Atlas of Peculiar 
Galaxies ^rp,|1966). It has been cla ssified as a SA (rs) a pec 



in RC2 (de Vaucouleurs et al. 1976| ), as a S O pec (Lynds & 



Furenlid 19 73|), as a n ellipticaTwith dust lanes (Hawarden et al. 



1981 i Bertolari987t| Yasuda et al.|1992l ), and as a BCG JThuan^ 
Martin 1981). Whit more et al.| ( |1990| ) considered III Zw 102 to 
be "related to polar-ring galaxies". Broad-band photometry was 



presented in|Cair6s et al.|(|2001a 



first studied by Demoulin 



b). Ill Zw 102 kinematics was 



(1969), which concludes, from opti- 



cal spectroscopy, that the galaxy rotates as a nearly solid body 
within the central 30''. Hi gas displays a more complicated be- 
havior, with the outer H i showing apparent counter-rotation with 
respect to the stellar and molecular gas components ( |Li et al.| 
[1953]). 

Results from 3D investigations of III Zw 102 have been 
published in Garcia-Lorenzo et al. (2008), where the central 
33.6" X 29.4" of the galaxy were mapped with a spa tial reso- 
lution of 2.7" spaxel"\ using the IFS integral, and in Moiseev 
2008| ), where the distribution and kinematics of the ionized gas 



was studied by means of Fabry-Perot scanning interferometry. 

Maps of III Zw 102 are shown in Figures |6] and [7] (col- 
umn 5). The VIRUS-P IFU 4'.' 2 fiber diameter corresponds to 
about 470 pc. 

In the continuum the galaxy displays outer circular 
isophotes, indicative of the presence of an underlying popula- 
tion of older stars, extending up to a radius of about 7 kpc. The 
central regions are more irregular, and have elliptical isophotes, 
with the major axis oriented in the northeast- southwest direc- 
tion. The presence of dust is clear seen on the continuum map: a 
dust lane crosses the galaxy in a southwest-northeast direction, 
almost parallel to the ellipses' major axis, and large dust patches 
are also evident at the south, almost parallel to the galaxy minor 
axis. In the narrow-band maps the morphology is more irregu- 
lar, with several filaments, bubbles and holes. The extended fila- 
ments, departing to the southeast and to the northwest, are very 
well traced in the [O ii] and HyS maps. 

The [O iii]/HyS ratio (see Figure [7]) is lower in the central part 
of the galaxy, with a minimum value o f ^ 0.22, in good agree- 
ment with the values published by Garcia-Lorenzo et al. I (120081). 
The excitation, like in Haro 1 , increases towar ds the g alaxy out- 
skirts, a behavior often seen in spiral galaxies ( Smith|1975| ). Just 



like Haro 1, III Zw 102 is a luminous object, which has been 
classified as a spiral by de Vaucoule urs et aLl ( |1976| ). 

The extinction map is highly inhomogeneous, and traces 
well the dust lanes detected in optical frames and color maps. 
The values are indicative of a considerable amount of dust. 

The stellar kinematics display an overall rotation around an 
axis oriented approximately northwest- southeast, with the north- 
east side approaching and the southwest receding, and a velocity 
amplitude of about 160 km s"^. The gas kinematics, which in 
this galaxy is better mapped by the HyS line (stronger than [O iii]) 
broadly agrees with the stellar kinematics in the inner region. 
Notable are the two features protruding from either side of the 
galaxy, and which, together with the knot 40" northeast of the 
galaxy center, seem to form a structurally and kinematically in- 
dependent system. These findings are in line with the results of 
the detailed analysis of this galaxy published by Moiseev ( 2008] ): 
III Zw 102 is made up of two distinct systems: an inner disk, 
where gas and stars co-rotate, and a polar ring orthogonal to it. 

The oxygen abundances, 12 -h log (0/H) = 8.44 in the inte- 
grated spectrum, and 12 -h log (0/H = 8.49 in the nuclear spec- 
trum, are much lower than those derived in > G arcia-Lorenzo et al.| 
( 2008) using the empirical calibrator N2 (De nicolo et al.|2002[ ), 
ranging from 8.81 to 8.84 in the circumnuclear regions and about 
8.85 for the nucleus. While the disagreement can be due partly to 
the diff'erent corrections for the underlying absorption at HyS, dis- 
crepancies of up to 0.4 dex between values obtaining using the 
P-method and N2 have already been reported ( |Shi et al.|2005 ). 



5. Summary 

We present results of an IFS analysis of five luminous BCG 
galaxies, based on VIRUS-P data covering the wavelength in- 
terval 3550-5850 A with a resolution of about 5 A FWHM. The 
mapped area is If 8 by If 8, each fiber corresponding to a circle 
of4'.'2onthe sky. 

For each galaxy, spaxel by spaxel, we modeled and sub- 
tracted the stellar continuum before fitting the most important 
emission-lines with a Gaussian. We produced an atlas of maps 
of the most relevant emission-lines in the spectra: [O ii] A3121, 
Hy, HjS, and [Oiii] /15007, of the continuum, of the extinction 
(Ay) and the excitation ([Oiii] ASOOl/HjS) ratios, as well as of 
the gas and the stellar kinematics. Integrated and nuclear spec- 
troscopic properties (line fluxes, extinction coefficient, and metal 
abundances) were also derived. From this study, we highlight the 
following results: 

- All galaxies show a regular morphology in the continuum, 
with outer elliptical or circular isophotes; the ionized gas 
displays a more complex pattern, with filamentary emis- 
sion and the presence of bubbles and holes. Remarkable are 
NGC 4670, where the gas emission shows two large exten- 
sions north and southwest of the center, and III Zw 102, 
where the two narrow features protruding almost symmet- 
rically from the main body, and a knot north of the galaxy 
center, form part of a polar ring. 

- The galaxies display diff'erent ionization patterns: in 

II Zw 33 and Mrk 314 the excitation ratio, [Oiii] ASOOl/Ufi 
traces the SF regions, as expected in objects ionized by 
stars. NGC 4670 has a somewhat complex pattern, with 
[O III] A5001/UJ3 peaking in the central [H ii] region and at 
the outer edge of the two large filaments. These three galax- 
ies have excitation ratios larger than one. In Haro 1 and 

III Zw 102, HyS is stronger than [Oiii], and the excitation 
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does not trace the SF regions, but has a minimum in the cen- 
tral galaxy regions (where the emission peaks) and increases 
towards the galaxy outskirts. 

- All galaxies possess significant amounts of dust, with inter- 
stellar reddening values up to Ay ~ 1.0. 

- A variety of kinematical behaviors are present in the sam- 
ple galaxies. II Zw 33 shows little or no rotation, while in 
Haro 1 gas and stars both display the same overall regular 
rotation pattern. The kinematics of NGC 4670 is typical of 
a barred galaxy, with S- shaped iso velocity contours in the 
gaseous component; stars co-rotate with the gas, albeit with 
a smaller gradient and amplitude. Mrk 314 has an irregular 
rotation pattern, with a drop in velocity in its southern lobe. 
The most complex object is III Zw 102, with two kinemati- 
cally independent systems: an inner disk, where gas and stars 
co-rotate, and a polar ring orthogonal to it. 

Acknowledgements. L. M. Cairos acknowledges the Alexander von Humboldt 
Foundation. VIRUS-P has been made possible by a generous donation from 
the Cynthia and George Mitchell Foundation. This paper includes data taken 
at The McDonald Observatory of the University of Texas at Austin. We thank 
J. N. Gonzalez-Perez and J. Falcon Barroso for stimulating discussions and 
tips on the data processing and analysis. N. Caon is grateful for the hospi- 
tality of the Leibniz-Institut fiir Astrophysik. This research has made use of 
the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet 
Propulsion Laboratory, Caltech, under contract with the National Aeronautics 
and Space Administration. This work has been funded by the Spanish Ministry 
of Science and Innovation (MCINN) under the collaboration "Estallidos" 
(grants HA2006-0032, AYA2007-67965-C03-01 and AYA2010-21887-C04-04). 
R Papaderos is supported by Ciencia 2008 Contract, funded by FCT/MCTES 
(Portugal) and POPH/FSE (EC). 



References 

Alter, L. H., ed. 1984, Astrophysics and Space Science Library, Vol. 112, Physics 

of thermal gaseous nebulae 
Arp, H. 1966, ApJS, 14, 1 

Bergvall, N. 2012, Star Forming Dwarf Galaxies, ed. P. Papaderos, S. Recchi, & 

G. Hensler, 175 
Bergvall, N. & Ostlin, G. 2002, A&A, 390, 891 

Bertola, F. 1987, in lAU Symposium, Vol. 127, Structure and Dynamics of 

Elliptical Galaxies, ed. R T. de Zeeuw, 135-144 
Blanc, G. A., Heiderman, A., Gebhardt, K., Evans, N. J., & Adams, J. 2009, ApJ, 

704, 842 

Burstein, D. & Heiles, C. 1982, AJ, 87, 1165 

Cairos, L. M., Caon, N., Garcia-Lorenzo, B., et al. 2007, ApJ, 669, 251 
Cairos, L. M., Caon, N., Garcia-Lorenzo, B., Vflchez, J. M., & Munoz-Tunon, 

C. 2002, ApJ, 577, 164 
Cairos, L. M., Caon, N., Papaderos, R, et al. 2009a, ApJ, 707, 1676 
Cairos, L. M., Caon, N., Papaderos, R, et al. 2003, ApJ, 593, 312 
Cairos, L. M., Caon, N., Vflchez, J. M., Gonzalez-Perez, J. N., & Munoz-Tunon, 

C. 2001a, ApJS, 136, 393 
Cairos, L. M., Caon, N., Zurita, C, et al. 2010, A&A, 520, A90 
Cairos, L. M., Caon, N., Zurita, C, et al. 2009b, A&A, 507, 1291 
Cairos, L. M., Vflchez, J. M., Gonzalez Perez, J. N., Iglesias-Paramo, J., & Caon, 

N. 2001b, ApJS, 133, 321 
Caon, N., Cairos, L. M., Aguerri, J. A. L., & Munoz-Tunon, C. 2005, ApJS, 157, 

218 

Cappellari, M. & Emsellem, E. 2004, PASP, 116, 138 

de Vaucouleurs, G, de Vaucouleurs, A., Corwin, Jr., H. G, et al. 1991, 
Third Reference Catalogue of Bright Galaxies, ed. de Vaucouleurs, G, de 
Vaucouleurs, A., Corwin, H. G, Jr., Buta, R. J., Paturel, G, & Fouque, P. 

de Vaucouleurs, G, de Vaucouleurs, A., & Corwin, J. R. 1976, in Second refer- 
ence catalogue of bright galaxies, 1976, Austin: University of Texas Press. 

Deeg, H., Duric, N., & Brinks, E. 1997, A&A, 323, 323 

Demoulin, M.-H. 1969, ApJ, 157, 69 

Denicolo, G, Terlevich, R., & Terlevich, E. 2002, MNRAS, 330, 69 
Dopita, M. A. & Sutherland, R. S. 1995, ApJ, 455, 468 
Doublier, V, Caulet, A., & Comte, G. 1999, A&AS, 138, 213 
Doubher, V, Comte, G, Petrosian, A., Surace, C, & Turatto, M. 1997, A&AS, 
124, 405 

FaneUi, M. N., O'Connell, R. W., & Thuan, T. X. 1988, ApJ, 334, 665 
Garcia-Lorenzo, B., Cairos, L. M., Caon, N., Monreal-Ibero, A., & Kehrig, C. 
2008, ApJ, 677, 201 



Gfl de Paz, A. & Madore, B. R 2005, ApJS, 156, 345 

Gfl de Paz, A., Madore, B. R, & Pevunova, O. 2003, ApJS, 147, 29 

Gordon, D. & Gottesman, S. T. 1981, AJ, 86, 161 

Grimes, J. R, Heckman, T., Aloisi, A., et al. 2009, ApJS, 181, 272 

Guseva, N. G, Papaderos, R, Izotov, Y. I., et al. 2003a, A&A, 407, 75 

Guseva, N. G, Papaderos, R, Izotov, Y. I., et al. 2003b, A&A, 407, 91 

Guseva, N. G, Papaderos, R, Izotov, Y. I., et al. 2003c, A&A, 407, 105 

Hawarden, T. G, Longmore, A. J., Tritton, S. B., Elson, R. A. W., & Corwin, Jr., 

H. G. 1981, MNRAS, 196, 747 
Hill, G. J., MacQueen, P. J., Adams, J., et al. 2007, in Bulletin of the American 

Astronomical Society, Vol. 39, American Astronomical Society Meeting 

Abstracts, 747 

Hill, G. J., MacQueen, R J., Smith, M. R, et al. 2008, in Presented at the 
Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 
7014, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 

Series 

Hill, G. J., MacQueen, R J., Tufts, J. R., et al. 2006, in Presented at the 
Society of Photo-Optical Instrumentation Engineers (SPIE) Conference, Vol. 
6269, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series 

Hunter, D. A. & Elmegreen, B. G. 2004, AJ, 128, 2170 
Hunter, D. A., van Woerden, H., & Gallagher, III, J. S. 1994, ApJS, 91, 79 
Izotov, Y. I., Guseva, N. G, Fricke, K. J., et al. 2010, A&A, 517, A90 
Izotov, Y. I., Schaerer, D., Blecha, A., et al. 2006a, A&A, 459, 71 
Izotov, Y. I., Stasihska, G, Meynet, G, Guseva, N. G, & Thuan, T. X. 2006b, 
A&A, 448, 955 

James, B. L., Tsamis, Y. G, Barlow, M. J., et al. 2009, MNRAS, 398, 2 

Klein, U., Wielebinski, R., & Thuan, T. X. 1984, A&A, 141, 241 

Kniazev, A. Y, Pustilnik, S. A., Grebel, E. K., Lee, H., & Pramskij, A. G. 2004, 

ApJS, 153, 429 
Kunth, D. & Ostlin, G. 2000, A&A Rev., 10, 1 

Li, J. G, Seaquist, E. R., Wrobel, J. M., Wang, Z., & Sage, L. J. 1993, ApJ, 413, 
150 

Loose, H.-H. & Thuan, T. X. 1986, in Star-forming Dwarf Galaxies and Related 

Objects, 73-88 
Lynds, B. T. & FurenHd, I. 1973, ApJ, 186, 445 

Marlowe, A. T., Meurer, G. R., Heckman, T. M., & Schommer, R. 1997, ApJS, 
112, 285 

Mas-Hesse, J. M. & Kunth, D. 1999, A&A, 349, 765 
Mas-Hesse, J. M., Kunth, D., Tenorio-Tagle, G, et al. 2003, ApJ, 598, 858 
Mazzarella, J. M. & Boroson, T. A. 1993, ApJS, 85, 27 
McCall, M. L., Rybski, R M., & Shields, G. A. 1985, ApJS, 57, 1 
Mendez, D. I., Cairos, L. M., Esteban, C, & Vflchez, J. M. 1999, AJ, 117, 1688 
Moiseev, A. V. 2008, Astrophysical Bulletin, 63, 201 
Moustakas, J. & Kennicutt, Jr., R. C. 2006, ApJS, 164, 81 
Murphy, J. D., Gebhardt, K., & Adams, J. J. 2011, ApJ, 729, 129 
Noeske, K. G, Papaderos, R, Cairos, L. M., & Fricke, K. J. 2005, A&A, 429, 
115 

Osterbrock, D. E. & Ferland, G. J. 2006, Astrophysics of gaseous nebulae and 
active galactic nuclei (Astrophysics of gaseous nebulae and active galactic nu- 
clei, 2nd. ed. by D.E. Osterbrock and G.J. Ferland. Sausalito, CA: University 
Science Books, 2006) 
Papaderos, R, Guseva, N. G, Izotov, Y. I., & Fricke, K. J. 2008, A&A, 491, 113 
Papaderos, R, Izotov, Y. I., Fricke, K. J., Thuan, T. X., & Guseva, N. G. 1998, 
A&A, 338, 43 

Papaderos, R, Loose, H.-H., Thuan, T. X., & Fricke, K. J. 1996, A&AS, 120, 
207 

Paturel, G, Petit, C, Prugniel, R, et al. 2003, A&A, 412, 45 
Peterson, C. J., Thonnard, N., Rubin, V. C, & Ford, Jr., W. K. 1978, ApJ, 219, 
31 

Pilyugin, L. S. 2000, A&A, 362, 325 
Pilyugin, L. S. 2001, A&A, 369, 594 
Pilyugin, L. S. & Thuan, T. X. 2005, ApJ, 631, 231 
Renka, L. L. & Cline, A. K. 1984, Rocky Mountain J. Math, 14, 223 
Sanchez-Blazquez, P., Peletier, R. R, Jimenez- Vicente, J., et al. 2006, MNRAS, 
371,703 

Schaerer, D., Contini, T., & Pindao, M. 1999, A&AS, 136, 35 
Schlafly, E. R & Finkbeiner, D. R 2011, ApJ, 737, 103 
Schlegel, D. J., Finkbeiner, D. R, & Davis, M. 1998, ApJ, 500, 525 
Schneider, S. E. & Salpeter, E. E. 1992, ApJ, 385, 32 
Searle, L. 1971, ApJ, 168, 327 

Shalyapina, L. V, Moiseev, A. V, Yakovleva, V. A., Hagen-Thorn, V. A., & 

Barsunova, O. Y. 2004, Astronomy Letters, 30, 583 
Shi, R, Kong, X., Li, C, & Cheng, R Z. 2005, A&A, 437, 849 
Smith, H. E. 1975, ApJ, 199, 591 

Stasihska, G. 2009, What can emission lines teU us?, ed. J. Cepa, 1 
Thuan, T. X. & Martin, G. E. 1981, ApJ, 247, 823 
Vacca, W. D. & Conti, R S. 1992, ApJ, 401, 543 



15 



L. M. Cairos et al.: Mapping Luminous Blue Compact Galaxies with VIRUS-P 

van Driel, W., Amaboldi, M., Combes, K, & Sparke, L. S. 2000, A&AS, 141, 
385 

Vanzi, L., Cresci, G., Telles, E., & Melnick, J. 2008, A&A, 486, 393 
Vorontsov-Velyaminov, B. A. 1959, in Atlas and catalog of interacting galaxies 
(1959), 

Walter, R, Brinks, E., Duric, N., & Klein, U. 1997, AJ, 113, 2031 
Whitmore, B. C, Lucas, R. A., McElroy, D. B., et al. 1990, AJ, 100, 1489 
Yasuda, N., Fujisawa, K., Sofue, Y., et al. 1992, PASJ, 44, 1 



16 



